Is there an attraction between spinons in the Haldane— Shastry model? 
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While the Bethe Ansatz solution of the Haldane-Shastry model appears to suggest that the 
spinons represent a free gas of half-fermions, Bernevig, Giuliano, and Laughlin (BGL) P, 0] have 
concluded recently that there is an attractive interaction between spinons. We argue that the dressed 
scattering matrix obtained with the asymptotic Bethe Ansatz is to be interpreted as the true and 
physical scattering matrix of the excitations, and hence, that the result by BGL is inconsistent with 
an earlier result by Essler We critically re-examine the analysis of BGL, and conclude that there 
is no interaction between spinons or spinons and holons in the Haldane-Shastry model. 
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The Haldane-Shastry model (HSM) HHIlllll plays 
a unique role among the integrable models of spin S — ^ 
chains. In what might be refered to as a brilliant the- 
oretical coup, Haldane and Shastry discovered indepen- 
dently in 1988 that a trial wave function proposed by 
Gutzwiller |^ in 1963 provides the exact ground state 
to a Heisenberg type spin Hamiltonian whose interaction 
strength falls off as the inverse square of the distance 
between two spins on the chain. If one imposes peri- 
odic boundary conditions (PBCs), and embeds the one- 
dimensional chain into a two-dimensional complex plane 
by mapping it onto the unit circle with the S ~ ^ spins 
located at complex positions rja = exp(i^Q!), where N 
denotes the number of sites and a — I, . . . , N , the Hamil- 
tonian 




possesses the exact ground state 
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^o{zi, . . . , zm) ^ Ylizi - Zjf Y[ Zj, (2) 

for N even, M = ^. The corresponding state vector is 
given by 

|vi'o)= Mzi, ■ ■■,zm) si- ■■ si, \U...i), 

{zi,...,Zm} 

where the sum extends over all possible ways to distribute 
the positions Zi of the up spins over the N sites. The 
model is fully integrable even for a finite number of sites; 
the algebra of the (infinite) number of conserved quanti- 
ties is generated by the total spin and rapidity operators 

N . N 

5 = 5:5., A=^^ Ih^^s^xS,) (4) 

^ 2 ^ 77a - 77/3 

which both commute with the Hamiltonian but do not 
commute mutually. The unique feature of the model. 



from a practical point of view, is that in addition to its 
amenability to solution by the asymptotic Bethe Ansatz 
(ABA) H m m 113, the ground state and many of 
the excited states (in principle, all the ones where the 
spins of the spinon excitations are fully polarized) can 
be written down in closed form, i.e., the wave functions 
are known explicitly. In particular, the wave function 
for an individual spinon excitation, which carries spin 
i but no charge, at site rja is constructed in complete 
analogy to the wave function for a quasihole excitation 
in a fractionally quantized Hall liquid : 
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^'a(zi, . . . , Zm) ^Y[iVa~ Zj) Y[{Z, - Zj)'^ Jl Zj (5) 
j=l i<j j=l 

where N odd, Af — The model may hence be used 

to illustrate the sense in which spinons are fractionally 
quantized excitations: the spin of the spinon is one-half, 
while the Hilbert space Q is built up from spin-flips, 
which carry spin one. 

On a more profound level, the model is unique in 
that there is no spin exchange between spinon excita- 
tions 0, which follows directly from the commutativity 
of A with i?HS- Furthermore, the spinon excitations of 
the model have or had been considered to constitute an 
ideal gas of half-fermions, that is, an ideal gas of parti- 
cles obeying fractional statistics [l4l Il5l | . This view has 
received strong support from Essler [3j, who calculated 
the dressed scattering matrix of the spinon excitations 
using the ABA, and found it to be S' = ±i. The fact 
that S does not depend on the spinon momenta implies 
that they are non-interacting or free; the phase i implies 
that they obey half-fermion statistics. This picture, and 
in particular the applicability of the ABA to the HSM, 
were commonly accepted until a few years ago. 

In 2001, this picture was challenged by Bernevig, Giu- 
liano, and Laughlin (BGL) "7, "^l, who investigated the 
nature of the spinon interaction by working out the wave 
functions for the spin-polarized two-spinon eigenstates 
explicitly. They found "clear evidence for a short-range, 
attractive interaction between spinons" 0|. Further- 
more, they "prove rigorously that this enhancement" — 
meaning a probability enhancement as the spinons are 



close together — "is responsible for the square-root sin- 
gularity in the dynamical spin susceptibility" 0, which 
has been evaluated exactly in the thermodynamic limit 
for the HSM by Haldane and Zirnbauer , and experi- 
mentally observed in KCuFa by Tcnnant et al. 0|. Ac- 
cording to BGL, "the experiments provide evidence that 
spinous do interact and that the spinon interaction is 
what determines the peculiar low-energy physics of spin- 
i antiferromagnetic chains" 0. BGL attribute the ap- 
parent contradiction between their results and the ABA 
result to the fact that "the interaction between spinous is 
encoded in the definition of the pscudomomenta" which 
label the Bethe Ansatz solutions In other words, they 
assert that it is a special feature of the ABA technique 
that the spinon excitations appear to be free, while there 
is in fact an attractive interaction between them. 

This line of reasoning may sound convincing at 
first sight. Indeed, in models like the Calogero- 
Sutherland |0, 0| or the Haldane-Shastry model, the 
long-range interaction of the particles or spins, respec- 
tively, is encoded in the definition of the pscudomomenta. 
The interacting degrees of freedom the Hilbert space is 
built up from, the particles or spin-flips, are mapped 
through a non-local and highly non-trivial transforma- 
tion into a new set of degrees of freedom, the pscudomo- 
menta, which do not interact. In a sense, in the HSM 
both the tail of the spin-flip terms S^S^ as well 
as the "potential energy" term S^Sp are encoded in the 
pscudomomenta. In the framework of the ABA, spinon 
excitations for the HSM correspond to fractional holes in 
the otherwise uniform distribution of pscudomomenta. 
Specifically, a pair of spinous is constructed by shifting 
the pscudomomenta quantum numbers li from integer to 
half-integer values or vice versa, and leaving the /^'s or 
pscudomomenta associated with the spinous unoccupied. 
The energy of the state is given by a sum of "kinetic en- 
ergies" of each occupied pseudomomentum, without an 
interaction between them. The information regarding the 
interaction between the original spins is no longer 
accessible in this framework. 

What is still accessible, however, is the information re- 
garding the energies of and the interaction between the 
spinon excitations. The energies of the spinous are given 
by the change in the kinetic energies associated with the 
occupied pscudomomenta as we shift them. The inter- 
action between the spinons is encoded in the way this 
shift in the pscudomomenta induced by one spinon is af- 
fected by the existence of another. In the spin one-half 
Heisenberg chain, for example, there is a rather compli- 
cated change or "screening" of the pscudomomenta due 
to an interaction between the spinons. In the HSM, by 
contrast, the creation of a spinon only induces a con- 
stant shift of the pscudomomenta, which implies that the 
spinons are free. The most reliable way to extract this 
information, however, is to calculate the spinon-spinon 
scattering matrix. If the ABA is applicable to the HSM 
at all, which is not garanteed a priori as the spin-spin 
interaction is long-ranged, the result by Essler quoted 



above unambigously confirms that the spinons are free. 

In the remainder of this Letter, we resolve the con- 
tradiction between the conclusions reached by Essler Q 
and BGL [l.,, I2|. The result is that we find no reason 
to doubt the applicability of the ABA, and completely 
agree with Essler's conclusions. We also agree with the 
explicit calculations of BGL, but do not agree with their 
interpretation of the calculations. In particular, their 
conclusion that there is an attraction between spinons or 
spinons and holons in the HSM, title to several publica- 
tions Ji, 1^, JJb , is incorrect. 

To begin with, it is worth noting that there is a physical 
reason to be suspicious of BGL's result. They conclude 
that there is an attractive interaction between spinons, 
but no bound state. If there was an arbitrarily weak 
attraction, however, it would presumably yield a bound 
state due to the Cooper instability j22j. Cooper's argu- 
ment was originally formulated for two electrons outside 
a completely occupied Fermi sphere, which are subject 
to an arbitrarily weak attraction. The argument is inde- 
pendent of the number of dimensions. The Fermi surface 
is only relevant in that it blocks certain states, and ren- 
ders the density of states available to the two electrons 
at the point where their kinetic energy is minimal (i.e., 
at the Fermi surface) finite. The Fermi statistics of the 
electrons accounts for the formation of a spin singlet, but 
is not essential to the instability; for example, one would 
also find a bound state if one were to use spinless bosons 
instead. The only subtlety involved in applying the ar- 
gument to spinons in the HSM is the half-fermi statistics 
of the spinons. It is not plausible to us, however, that 
this statistical interaction would preclude the pairing, as 
there is not even an angular momentum barrier associ- 
ated with the statistical interaction in one dimension. 

Let us now critically re-examine the arguements pre- 
sented by BGL. We begin with a review of their analysis, 
and then explain why we disagree. 

BGL construct exact two-spinon eigenstates for the 
HSM starting from basis states with the two spinons lo- 
calized at sites rja and ?7^, 

M M M 

J = l i<j J = l 

(6) 

where M — ^ ^ ^ denotes the number of up or down spins 
condensed in the uniform singlet sea. The momentum 
space basis states are obtained by Fourier transformation. 



^mn{Zi, Zm) = 2^ — *Q/3(Zl, . . . , Zm), 

a, (3 

(7) 

where M > m > n > 0. For m or n outside this range, 
^mn will vanish identically, refiecting the overcomplete- 
ness of the position space basis Acting with the 
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Haldane-Shastry Hamiltonian on Q yields 



^HS l^mn) — Ej, 



mn |,T, \ 



1=1 



where 



Emn = -J^^^- 



19 24 



J f2n 



m| — — 1 — m 



N 



N 
l~n 



(8) 



(9) 



Im =-aim{M-m,n), and Vi""' ^ - ^ {m - n + 21) . 

Since the "scattering" of iJns acting on the non-ortho- 
gonal basis states l^'mn) only occurs in one direction, 
increasing the difference m — n while keeping the "total 
momentum" m + n fixed, the (unnormalized) eigenstates 
of iJns have energy eigenvalues Emn and are of the form 



singularity at the lower threshold frequency of the two- 
spinon continuum, which is a characteristic feature of 
spin-i chains [2^ . 

As already mentioned, we completely agree with these 
calculations. We disagree, however, with BGL's interpre- 
tation of the results as evidence for a spinon attraction. 

The first argument given by BGL in favor of a spinon 
interaction is based on a plot of |pmn(e*^)P for m — Af, 
71 = 0, as a function of 9. They interpret |pmn(e**)P as 
probability for finding the spinous at a distance 9 along 
the circle from each other, and show it to be strongly 
enhanced at small 9. The problem with the argument 
is that, as one can easily see from (fTT)) . the PmniVa-^ys 
are the coefficients in the expansion of the overcomplete 
basis states l^'a/s) at fixed a, (3 in terms of |<i>m„). Due 
to this overcompleteness, the PmniVa-pYs as functions 
of rja-p have no direct physical interpretation. The ac- 
tual relative spinon-spinon wave function (fmniila-p) for 
given m and n provides the coefficients in 



(10) 



1=0 



with a™" = 1. A recursion relation for the coefficients 
a™" is easily obtained from (jHJ. Combining H1U|I and ||7J), 
one thus obtains an expansion of the exact energy eigen- 
states |$m„) in terms of localized spinon states I^'q/s)- 

In a technically truly remarkable analysis, BGL have 
further succeeded in obtaining the coefficients Pmnixia-p) 
in the inverse expansion 

M m 

l*"/3) = E E(-1)'"^"C^^ Pran{llo.-p) \<^mn) (H) 
m— n— 

of the localized spinon states in terms of the energy 
eigenstates by solving a hypergeometric differential equa- 
tion. Since a spin-flip S~ acting on the Haldane-Shastry 
ground state yields a state with a pair of spinous local- 
ized at 77a, S~ l^'o) = rja {^aa) , the expansion of S"^ |^'o) 
in terms of is determined by Prjm(l), 
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M m 
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OPmn(l) l^'mn) , (12) 



where q = and the Kronecker-(5 is defined modulo 
N. The explicit expression for Pm„(l) enabled BGL to 
calculate the dynamical spin susceptibility (DSS) 



X,(c^) = -Im(*ol^^ 



1 



uj~iHiis-Eo)+iO 



.(13) 

for finite chains as well as in the thermodynamic limit, 
thus providing an alternative derivation of the Haldane- 
Zirnbauer formula The DSS shows a square-root 
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|$mn) = E E 'PmniVa-p) ' (??a+/3)" 
a=l 13=1 



i+n |\J> 



(14) 

It is easily seen from ()10|l and Q that a possible choice 

for fmnilla-f}) is 

VmniVc-p) = E «™ • iVc-p)^^ mc.p)\\. (15) 
/=0 

Depending on m and n, one finds that <Pmn(e*^) is some- 
times enhanced and sometimes suppressed for small 9, 
but even if there was a clear enhancement, it would not 
allow for a conclusion regarding a spinon attraction. The 
reason is simply that the overcompleteness of the ba- 
sis states l^'a/s) implies that ipmniVa-is) is not uniquely 
determined, i.e., there are infinitely many choices for 
VmniVa-p) which yield the same \^rnn) in C3- 

The second argument of BGL is that the last term 
in the energy © of the two-spinon state \^mn) repre- 
sents "a negative interaction contribution that becomes 
negligibly small in the thermodynamic limit" jlj. The 
problem here is that BGL identify the momenta qm and 
g„ of the individual spinous according to 

■n 2tt f l\ 7r27r/ l^,^ 

9™=o-T7 "^+oh 9« = 77-l^ "+o . (16) 
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and interpret the two preceding terms in ^ as the kinetic 
energies of the individual spinons. The correct identifi- 
cation of the spinon momenta for m > n, however, is 



TT 2tt / 3 

2- AT 



TT 27r/ 1\ , , 

2-Ar^^+4)' ^''^ 



which implies that the kinetic energy of the spinons is 
given by all three terms in the square bracket in 
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With E{q) = I [(§ )^ - g2] , one finds 

Emn ^ ~J^(^N+^~^^ + E{q^) + E{q,,). (18) 

The alleged spinon interaction term has disappeared. 
Physically, the shift between qm and qn by one-half of 
a momentum spacing ^ is nothing but a manifestation 
of the half-fermi statistics of the spinous. While the al- 
lowed values for the total momenta g„i + Qn are those for 
PBCs, the allowed values for the difference in the mo- 
menta qm — Qn are those for anti-PBCs, i.e., PBCs with 
the ring threaded by a flux tt. 

Finally, BGL claim to prove that the enhancement of 
|Pmn(e**)P they find when plotting it as a function of 
the spinon separation 9 is responsible for the square- 
root singularity in the DSS 0. Their proof then con- 
sists of the derivation of the Haldane-Zirnbauer formula 
sketched above. In their longer paper 2], BGL conclude 
that their "analysis definitely proves that the square-root 
sharp edge on top of the broad spectrum is nothing but 
the interaction between spinous" , and say that "this re- 
sult is of the utmost importance, since it represents a way 
to experimentally test the interaction among spinous in 
one dimension" . 

There are several problems attached to this line of rea- 
soning. First, the coefficients |pmn(e'^)P cannot be inter- 
preted as a probability as a function of the spinon separa- 
tion 9, as explained above. Second, it is not Pmnis^^) as 
a function of 9 for fixed m and n which enters the deriva- 
tion of the Haldane-Zirnbauer formula, but Pmn(l) as a 
function of m and n, as one can directly see from (|12|l . 

The square-root singularity in the DSS is not due to an 
alleged spinon attraction, but a general consequence of 
the fractional quantization of spin excitations in spin-i 
chains. The position space basis for these fractional ex- 
citations, the spinons, is necessarily overcomplete. The 
local creation of two spinons through a spin-flip is not 
equivalent to a creation of all two-spinon energy eigen- 
states with the same relative weight. The process rather 
creates predominantly spinons with lower energies, which 
is reflected in the square-root singularity in the DSS. 
Since the fractional quantization of spin excitations is 
a generic feature of spin-i chains, the square-root singu- 



larity in the DSS is generic as well. It exists in the HSM, 
where spinons are free, but also in the Heisenberg model, 
were spinons are interacting. With the experimental ob- 
servation of the square-root singularity in KCuFs, Ten- 
nant et al. 17] have observed fractional quantization in 
spin chains. They did not observe a spinon-spinon at- 
traction. 

In the context of this analysis, it is worthwhile to men- 
tion a curiosity of the two-spinon eigenstates. In the 
evaluation reviewed above, BGL obtained the coefficients 
a™" by explicitly solving the Sutherland equation © us- 
ing the Ansatz (|10() . The coefficients a™" were hence 
determined by the Hamiltonian, and appear to contain 
information inflicted on the system by the Hamiltonian. 
In principle, this could include information regarding an 
interaction between spinons. 

In fact, however, the Hamiltonian is not even required 
in determining the coefficients a™". If we wish to con- 
struct an orthogonal basis \^mn) according to (|10() with 
a™" = 1 from the non-orthogonal basis |^mri)i the over- 
laps (^'mnj^'m'n') for all TO, Tz, m', n' completely deter- 
mine all the coefficients aj"", as the reader will be able 
to verify easily for himself. The coefficients a™" as well as 
PmniVa), and therefore also the "scattering amplitudes" 
V™" in ||SJ), hence contain no information except the one 
regarding the Hilbert space structure of the fractionally 
quantized excitations. Accordingly, it seems impossible 
as a matter of principle to reach a conclusion regarding 
an interaction between the spinons by studying these co- 
efficients. 

In conclusion, we have shown that the spinons in the 
HSM represent an ideal gas of half-fermions, and thereby 
dispersed all evidence that the ABA might not be ap- 
plicable to the model. An analysis similar to the one 
presented here shows that there is likewise no interaction 
between spinons and holons in the HSM. The conclusions 
drawn by BGL with regard to this question [l9l f20| are 
likewise incorrect. 
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